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ABSTRACT 

We investigate in detail the hypothesis that low surface brightness galaxies (LSB) differ 
from ordinary galaxies simply because they form in halos with large spin parameters. 
We compute star formation rates using the Schmidt law, assuming the same gas infall 
dependence on surface density as used in models of the Milky Way. We build stellar 
population models, predicting colours, spectra, and chemical abundances. We compare 
our predictions with observed values of metallicity and colours for LSB galaxies and 
find excellent agreement with all observables. In particular, integrated colours, colour 
gradients, surface brightness and metallicity match very well to the observed values 
of LSBs for models with ages larger than 7 Gyr and high values (A > 0.05) for the 
spin parameter of the halos. We also compute the global star formation rate (SFR) 
in the Universe due to LSBs and show that it has a flatter evolution with redshift 
than the corresponding SFR for normal discs. We furthermore compare the evolution 
in redshift of [Zn/H] for our models to those observed in Damped Lyman a systems 
by Pettini et al. (1997) and show that Damped Lyman a systems abundances are 
consistent with the predicted abundances at different radii for LSBs. Finally, we show 
how the required late redshift of collapse of the halo may constrain the power spectrum 
of fluctuations. 

Key words: galaxies: formation - galaxies: evolution - galaxies: spiral - galaxies: 
stellar content. 
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1 INTRODUCTION 

Late type low surface brightness galaxies are found to be 
a significant fraction of the total galaxy population in the 
Virgo cluster (Impey, Bothun & Malin 1988), in the For- 
nax cluster (Irwin et al. 1990), and in the field (McGaugh, 
Bothun & Schombert 1995; Sprayberry, Impey & Irwin 1996; 
Sprayberry et al. 1997). A theory of galaxy formation should 
therefore account for the existence of LSBs. On the other 
hand, LSBs are particularly useful because they are sim- 
pler objects than high surface brightness (HSB) galaxies: 
they have relatively little present day star formation and 
little dust (reddening). Therefore it is easier to model their 
stellar population and star formation history. They are also 
more dark matter dominated than HSBs, and therefore par- 
ticularly suitable for probing the structure of dark matter 
halos. In order to compare theories of galaxy formation with 
observational surveys of galaxies it is important not only to 
quantify observationally the abundance of LSBs, but also to 
understand their nature. For instance, if LSBs are made of 
very young stellar populations, as often claimed in the litera- 
ture, they may be irrelevant in the picture of the Universe at 



high redshift, while they could be a substantial component 
of that picture if their stellar populations are old. 

The existence of LSB galaxies can be readily understood 
if it is assumed, following Fall & Efstathiou (1980), that the 
specific angular momentum of the baryons is approximately 
conserved during their dissipation into a rotationally sup- 
ported disk, and that the disk length-scale is therefore re- 
lated to the angular momentum of the dark matter halo. 
The low surface brightness is the consequence of the low 
surface density of the disk, which is due to the larger spin 
parameter of the dark matter halos of LSBs relative to the 
spin parameter of the halos of HSB disks (e.g. Dalcanton, 
Spergel & Summers 1997; Mo, Mao & White 1997; Jimenez 
et al. 1997). 

However, while the surface density of a disk can be sim- 
ply related to the mass and spin parameter of its dark matter 
halo, its surface brightness, and therefore its mass to light 
ratio (M/L), depends on the stellar population and could in 
principle vary with galaxy type and age. For instance, M/L 
must have a significant time dependence due to the history 
of star formation in the disk and to the continuous death of 
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Figure 1. Initial surface density for an isothermal sphere and 3 
different masses. The solid lines correspond to A = 0.1 to 0.01 
(from top to bottom) in steps of 0.01. The circular velocity of the 
halo (dashed line) is also plotted. 



massive stars. Moreover, M/L in a given photometric band 
changes with time also as a consequence of the colour evolu- 
tion of the stellar population, which is sensitive also to the 
history of the chemical enrichment of the disk. Therefore, a 
comparison between a simple model of disk formation and 
the observations cannot be made without an appropriate 
model for the stellar population in the disk that includes 
self-consistently the chemical evolution of the population. 

In this paper we model self-consistently the chemical 
and photometric evolution of the stellar population of LSB 
disks, formed in dark matter halos described by isothermal 
spheres. The aim of this paper is twofold: 

(i) show that only one parameter, namely the spin param- 
eter A of the dark halo, can explain the surface brightness 
of LSBs. 

(ii) use LSB disk galaxies to study the formation and evo- 
lution of all disk galaxies. 

The second point is motivated by the first, that is by 



the fact that LSB disks are in fact so similar to HSB disks. 
On the other hand, LSBs are more dark matter dominated 
than HSBs, and therefore better described by a very simple 
model where only the gravitational field of an isothermal 
halo is considered. 

Our assumptions are as follows: 

(i) The specific angular momentum of baryonic matter is 
the same as for dark matter 

(ii) Gas settles until centrifugally supported, in a given 
dark matter halo 

(iii) The star formation rate is given by the Schmidt law 

(iv) The gas infall rate is assumed to be the same function 
of total surface density as used in models of the Milky Way 

We first use the density profile of the halo to compute 
the surface density of the settling disk. We then use the 
Schmidt formation law and an infall rate that reproduces the 
observed properties of the Galaxy in conjunction with the 
chemical evolution models by Matteucci & Francois (1989) 
to compute the star formation rate at several radii of the 
disk and the evolution of several chemical species (H, D, 
He, C, N, O, Ne, Mg, Si, S, Ca, Fe and Zn). We then pro- 
ceed to compute spectra, integrated colours, colour profiles 
and surface brightness for two different values of the spin 
parameter of the halo. 

We first present in the next section the simple non-self- 
gravitational disk model and the non-singular halo model. 
We describe in section 3 the set of synthetic stellar popula- 
tion models used to predict the spectra and colours of the 
stellar population. In section 4 we model the star forma- 
tion in the disk with a Schmidt law and a gas infall law, 
assuming the same model of the chemical evolution of the 
Galaxy as Matteucci & Francois (1989). Finally, in section 
5 we compute synthetic spectra and photometric properties 
of the stellar population in the disk, using the stellar evo- 
lution tracks of JMSTAR15, stellar atmosphere models by 
Kurucz (1992) and j0rgensen (private communication) and 
the chemical evolution models built for LSBs (see section 4). 
We discuss the results in section 6. 

The most important results of this work are: 

(i) Observed colour profiles, chemical abundances and 
surface brightness profiles for LSBs are well fitted if they 
are assumed to have a spin parameter for its halo higher 
than HSBs. 

(ii) LSBs are not young objects, as often claimed in the 
literature, since their colours are well fitted by old (> 7 Gyr) 
stellar populations. 

(iii) There is discrepancy between the photometric age of 
the galaxies and the age of formation of their halo, which in- 
dicates that the star formation can start about 2 Gyr before 
the halo is formed. This discrepancy is reduced to 1 Gyr if 
the Hubble constant is assumed to be Ho = 65 kms _1 Mpc _1 
and completely removed if the Universe is open (fi < 0.3) 
or has a significant vacuum energy contribution (A > 0.6). 

(iv) The earliest stellar populations of LSBs could be 
present in the high redshift universe in a similar propor- 
tion to HSBs as they are now. Their colours at high redshift 
(as the colours of all disks) are about 1 mag bluer than at 
low redshift. 

In a previous paper (Padoan, Jimenez & Antonuccio- 
Delogu 1997) we have shown that LSBs are not necessar- 
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Figure 2. Contour plot of the value of the central surface density (in Mqpc 2 ) for different values of the halo mass and spin parameter. 
Left panel: halo formation redshift z = 0.7. Right panel: halo formation redshift z = 2.1. 



ily young and un-evolved systems (e.g. McGaugh & Bothun 
1994; De Blok, Van Der Hulst & Bothun 1995). We used the 
simplest possible model, a burst of star formation, to deter- 
mine a lower limit to the age of LSBs studied by De Blok, 
Van Der Hulst & Bothun (1995). In the present work we 
considerably improve our previous model because we use a 
more realistic continuous star formation process, we include 
self-consistently the chemical evolution, and we connect the 
disk model to the cosmological scenario using the spheri- 
cal collapse model (Gunn & Gott 1972). This more detailed 
model confirms our previous result that the blue LSB disks 
in the sample by De Blok, Van Der Hulst & Bothun (1995) 
are not un-evolved objects collapsed at late times from low 
initial over-densities (McGaugh & Bothun 1994; Mo, Mc- 
Gaugh & Bothun 1994), but rather normally evolved disk 
galaxies. 

Throughout the paper the present day value of the Hub- 
ble constant is assumed to be Ho = 75 kms _1 Mpc~ 1 . 



2 THE HALO AND DISK MODELS 

In this work we model LSB galaxies as non-self-gravitating 
disks inside isothermal dark matter halos, following Mo, Mao 
& White (1997). Isothermal spheres are found to present a 
rough description of halos in numerical simulations (Frenk 
et al. 1985; Frenk et al. 1988; Quinn, Salmon & Zurek 1986; 
Navarro, Frenk & White 1997), and have been previously 
used in galaxy formation models by Kauffmann, White & 
Guiderdoni (1993) and Cole et al. (1994). 

We briefly summarize the non-self-gravitating disk 
model as presented in Mo, Mao & White (1997). In the 
spherical collapse model (Gunn & Gott 1972), the halo mass 
M and its circular velocity V c , at redshift z, are: 



M = 



10GH(z) 



(1) 



where H(z) is the Hubble constant at redshift z: 

H(z) = H [n A>0 + (i-n A , ~Qo)(i + z) 2 + n i) (i + z) 3 ] 1/2 (2) 



where fio, fiA.o and Ho are the values at z = 0. Note that 
the halo is characterised by a circular velocity V c , but there 
is no implication that the halo is rotating. V c is the speed 
required for centrifugal support. In general, the halo is not 
expected to be rotationally supported. 

The disk is assumed to be thin and in centrifugal bal- 
ance, to have an exponential surface density profile, and an 
angular momentum and a mass that are a fraction /j of the 
angular momentum and mass of the dark matter halo. The 
value of /t, is here assumed to be equal to the ratio of the 
baryonic to total matter density: 



fb — fi^o/fio 



(3) 



The disk mass Ma, surface density Ed, and scale-length Rd 
arc: 



M A 



27rE i?d 



WGH(z 

E d (R) = E ex.p(-R/R d ) 

R _ XVc 
d 10V2H{z) 

where Eo is the central surface density 

Wf b V c H(z) 
E ° = -kGX 2 

and A is the spin parameter of the halo: 
A = J|£| 1/2 G _1 M- 5/2 



(4) 
(5) 
(6) 

(7) 
(8) 



and E is the total energy of the halo. 

In Fig. |l] surface density profiles are plotted for z — 
0.7 for three different values of the halo mass. In each plot 
the profiles for different values of A are shown, with 0.01 < 
A < 0.1. Also shown is the circular velocity. Models with 
larger A have lower central surface density and larger disk 
scale-length. Fig. ^| is a contour plot of the value of the 
central surface density for different values of the halo mass 
and A. As can be seen in equation (Q), Eo oc M 1//3 A -2 . A 
variation in A between 0.04 and 0.1 corresponds therefore to 
a change in Eo of a factor 6.25, which is equivalent to about 
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2 magnitudes in surface brightness, from the Freeman law 
value of /is(0) = 21.6 mag arcsec -2 (Freeman 1970), to 
jtts(O) = 23.6 mag arcsec -2 , if M/L of the baryons does 
not depend on A. On the other hand, the disk scale-length 
scales like Rd oc A, so a change in Rd of a factor 2.5 can be 
achieved. Since the halo is isothermal (p oc R~ 2 ), its mass 
grows proportionally to the distance from the center, and 
the ratio of dynamic mass, Md yn , to disk mass, is therefore 
proportional to the extension of the disk, 

M dyn /M dtsk oc A (9) 

because more extended disks contain a larger fraction of the 
dark matter halo than less extended disks, and the dark halo 
dominates the gravitational potential. This explains the fact 
that LSBs typically have M dyn /L about twice as large as 
HSBs, while they still obey the same Tully-Fisher relation as 
HSBs (Zwaan et al. 1995), as long as the M/L of the baryons 
is not significantly different in the two classes of galaxies: 
equations (1) and (6) imply V c 3 oc Md yn /X oc M d i S k- 



2.1 Isothermal Halos and Disk Properties 

In this work we refer to the sample of LSB disk galax- 
ies by De Blok, Van Der Hulst & Bothun (1995), for the 
photometry, and by De Blok, McGaugh & Van Der Hulst 
(1996), for the HI profiles and rotation curves. The HI 
data for U0128, U1230, U5005, U5750, and U6614 are from 
Van Der Hulst et al. (1993). Photometric data for the 
galaxies F563-V2 and F571-8 are taken from McGaugh & 
Bothun (1994) and De Blok & McGaugh (1997) respec- 
tively. Only 14 galaxies with HI rotation curves that clearly 
flatten have been used for the discussion in this section. 
The mean central surface brightness of the 14 galaxies is 
A's(O) = 23.5 mag arcsec -2 , the mean disk scale-length 
Rd — 4.7 kpc (Ho = 75 kms -1 Mpc -1 ), and the mean circu- 
lar velocity V c = 104 kms -1 . 

In Fig. ^ (upper panel) the HI disk scale-length of the 
models, Rd, is plotted versus the circular velocity, for dif- 
ferent values of formation redshifts z and spin parameter A. 
Also plotted are the 14 LSB galaxies (squared symbols) and 
two theoretical models for which the chemical and spectro- 
photometric evolutions are calculated in this work (aster- 
isks). The scale-length of the galaxies is the one measured 
in the B photometric band. The plot shows that, for a range 
of A that should be appropriate for LSBs, the halos of the 
galaxies in the sample should turn-around in the redshift 
range 0.1 < z < 1.2. 

A different way to estimate the size of the galaxies in 
this sample is to use the radius, = Ri, at which the HI sur- 
face density is equal to 1 Mqpc -2 . The value of Ri is plotted 
versus the circular velocity in Fig. ^ (middle panel) for the 
theoretical models and for the galaxies in the sample. The 
two models discussed in this work are also marked (aster- 
isks). This plot confirms the result of the previous one: the 
halos of the LSB galaxies in this sample are formed at a 
redshift in the range 0.1 < z < 1.2. Notice that the radius 
Ri is typically very large compared to the optical disk of 
the galaxy, so that virtually no star formation has occurred 
at that radius. This is also suggested by the fact that the 
HI profiles look very close to exponential around Ri, while 
most of the gas has clearly been lost into stars at inner 
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Figure 3. Upper panel: The HI disk scale— length of the models, 
Rd, versus the circular velocity, for different values of formation 
redshifts z and spin parameter A. The squared symbols are 14 LSB 
galaxies, and the asterisks the two theoretical models for which 
the chemical and spectro-photometric evolutions are calculated 
in the paper. The scale— length of the LSB galaxies is the one 
measured in the B photometric band. The plot shows that, for a 
range of A that should be appropriate for LSBs, the halos of the 
galaxies in the sample should turn-around in the redshift range 
0.1 < z < 1.2. 

Middle panel: The radius, = Ri , at which the HI surface density is 
equal to 1 Mqpc -2 , versus the circular velocity, for the theoretical 
models and for the galaxies in the sample (squared symbols). The 
asterisks are again the two models discussed in the paper. 
Bottom panel: The disk scale— length is plotted versus Hi. All 
symbols are like in the previous two panels. The galaxy that ap- 
pears inconsistent with the model is F583-1. 

radii. Therefore, the value of Ri is suitable for a direct com- 
parison between the model prediction and the observations. 
Conversely, the HI profiles are flattish or with holes at inner 
radii, which suggests that most of the gas has turned into 
stars inside Rd and therefore the value of Rd can also be 
compared directly with model predictions. 

The strong correlation between the value of the B band 
disk scale-length and the value of Ri, relative to the mod- 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy Formation and Evolution: Low Surface Brightness Galaxies 5 





Vc/ikms" 1 ) 


M/(l.elOM ) 


A 


2 


F561-1 


52 


3.6 


0.09 


0.1 


F563-1 


111 


26.4 


0.07 


0,1 


F563-V1 


30 


1.2 


0.07 


0.0 


F563-V2 


115 


20.0 


0.05 


0.8 


F568-1 


119 


7.9 


0.30 


2.5 


F568-V1 


124 


16.0 


0.10 


1.4 


F571-8 


133 


26.3 


0.11 


1.0 


P^71 VI 


t o 


i i 
i.i 


n 19 

U.lo 


1 n 

1 .u 


(F583-1 


85 


300.0 


0.001 


-0.8) 


U0128 


131 


26.6 


0.14 


0.9 


U1230 


102 


39.0 


0.05 


0.0 


U5005 


99 


22.9 


0.06 


0.2 


U5750 


75 


7.5 


0.14 


0.1 


U6614 


200 


178.5 


0.09 


0.2 


mean: 


105 


29.2 


0.11 


0.7 



Table 1. Circular velocity, halo mass, spin parameter and halo 
formation redshift from equations (10)-(12), for 14 LSB galaxies 
with rotation curves that clearly flatten. F583-1 has no solution 
and is not included to compute the mean values. 
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Figure 4. The absolute blue magnitude of the galaxies is plotted 
versus the estimated halo mass (from Table 1) and the value of 
V^Rd, that should be proportional to the halo mass, contained 
inside the radius R d . In both plots the absolute magnitude grows 
with the mass estimator. 
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Figure 5. Upper panel: The central surface brightness is plotted 
versus A. 

Middle panel: The central surface brightness is plotted versus the 
central surface density computed consistently with M, A, and z 
of each galaxy. 

Bottom panel: The central surface brightness is plotted versus 
V c 2 /R d . 



els, is also shown in the bottom panel of Fig. where Rd is 
plotted versus Ri. All galaxies but one are again described 
by formation redshift in the range 0.1 < z < 1.2. The galaxy 
that appears inconsistent is F583-1. This galaxy has an ex- 
tremely small scale-length Rd = 1.6 kpc (De Blok & Mc- 
Gaugh 1997), given its circular velocity V c = 85 kms - . 

A non-self-gravitating disk model is characterized by 
three parameters: the halo mass M, the spin parameter A, 
and the formation redshift z. For each galaxy in the sam- 
ple three parameters are known: the circular velocity V c , 
the disk scale-length Rd, and Ri. Using these three mea- 
sured quantities, the unknown parameters M, A and z can 
be determined for each galaxy. Using the equations above 
one gets: 



M 



h 



exp(Ri/R d ) 



(10) 
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Figure 6. Rotation curves and density profiles for the JHHP model. The solid lines correspond to different values for A (0.07 to 0.03, 
from top to bottom) and the dashed line is the corresponding critical surface density to form stars according to the Toomre/Kennicutt 
criterion. The predicted initial surface density for the disk is very similar to the isothermal sphere model. 



A = 



ftV c 2 



V2tySGR c 



exp(-i?i/i? d 



f b V?exp(-Ri/R d 



20irSGR%Ho 



2/3 



(11) 

(12) 



galaxies. The observed parameters of F583-1 are apparently 
inconsistent with the model, as mentioned above. 

In Fig. ^| the absolute blue magnitude of the galaxies 
is plotted versus the estimated halo mass (from Table 1) 
and the value of V^Rd, that should be proportional to the 
halo mass, contained inside the radius Rd- In both plots the 
absolute magnitude grows with the mass estimator. 

In Fig. U the central surface brightness is plotted ver- 
sus A (upper panel). There is no clear trend in the plot, 
as could be expected since the central surface density, and 
therefore the central surface brightness of a disk depends 
on both the spin parameter and the mass. When the cen- 
tral surface density of the disk is appropriately computed 
(equation (7)) consistently with the values of M, A, and z 
estimated for each galaxy (Table 1), a clear correlation be- 
tween the central surface brightness and the central surface 
density can be seen, as shown in Fig. ^ (middle panel). The 



observed central surface brightness grows approximately lin- 
early with the estimated central surface density, £d(0), for 
Sd(0) > 50 Mope -2 . If the halo formation redshifts and 
the spin parameters were unknown, the mass could only be 
estimated as something proportional to V^R d , and the cen- 
tral surface density as proportional to /Rd- In the bottom 
panel of Fig. |j the central surface brightness is plotted ver- 
sus V?/R d . There is no clear trend in the plot, which shows 
that a correct estimate of mass and surface density requires 
a knowledge of A and z for each galaxy, and not only of 
its circular velocity and size. The middle plot of Fig. M is 
therefore a strong indication that i) the surface brightness 
correlates with the surface density, and ii) the surface bright- 
ness is determined by the spin parameter, for any given halo 



2.2 Non-singular halo 

In Jimenez et al. (1997) a model (JHHP model hereafter) 
was developed for computing the final settling radius and 
surface density of the gas of a rotationally supported disk. 
In this model it was assumed that the gas component has 
the same specific angular momentum as the dark matter, 
and the dark matter profile found in simulations by Navarro, 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy Formation and Evolution: Low Surface Brightness Galaxies 7 






Figure 7. The model predictions for the JHHP model are compared with a selected sample of LSBs from De Blok, Van Der Hulst & 
Bothun (1995). The rotation curves are not well fitted in the inner parts of the halo. The right three panels show the observed HI surface 
density (De Blok, McGaugh &; Van Der Hulst 1996) compared to the previously computed E c using the Toomre/Kennicutt instability 
criterion. The agreement between both curves is good beyond 2 Kpc. 



Frenk & White (1997) was used. A detailed description of 
the model can be found in Jimenez et al. (1997). Using the 
above we have computed several disk models for different 
values of the mass and spin of the halo. Fig. ^ shows the 
predicted rotation curves and initial surface densities for 3 
different masses and 5 values for A: 0.03, 0.04, 0.05, 0.06 and 
0.07 (from top to bottom). Also plotted is the critical surface 
density (dashed line) to form stars according to the Toomre 
(1964) stability criterion (see also Kennicutt (1989)). If the 
surface density of the disk is above critical, then it will form 
stars, otherwise it will not. 

In Fig. we compare the predictions of the JHHP disk 
model with observed rotation curves and HI surface den- 
sities for F568-1, F568-V1, F563-1, F571-V1 and F563-V1. 
The central region of the rotation curves in all cases does not 
agree well with the model predictions since the observed cen- 
tral rotation curves are shallower than the predicted ones. 
The three right panels of Fig. |^ compare the observed HI 
surface density with the predicted critical surface density. In 
all cases the observed HI surface density and the predicted 
critical density agree reasonably well beyond 3-4 kpc show- 
ing the validity of the Toomre criterion, while in the central 
regions the observed HI surface density is much lower than 



the predicted critical density, pointing towards global insta- 
bilities controlling star formation in the central regions. 

The important point to notice here is that the JHHP 
disk model and the isothermal sphere predict roughly the 
same initial surface density for a rotationally-supported disk 
(as can be seen by comparing the bottom-right panel from 
Fig. |^ and the middle panel from Fig. ^). This shows that 
the initial surface density is a robust prediction and since, 
as it will be shown in the next two sections, it controls what 
the colours (and therefore ages) of an LSB will be, we regard 
our ages and colours as robust predictions. 

In order to compare a disk model with the photomet- 
ric properties of LSB galaxies, it is necessary to calculate 
an appropriate stellar population model, based on the pre- 
dicted surface density. It is not possible to calculate models 
of galaxies without a very detailed and up-to-date treat- 
ment of stellar evolution and atmosphere models, for the 
simple reason that galaxies are made of stars. The stellar 
population model must be based on the computation of the 
chemical evolution of the galaxy, consistent with a reason- 
able star formation history, because the colours of a stel- 
lar population depend strongly on its chemical composition, 
and the chemical composition of the star forming gas is de- 
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Figure 8. The Zero Age Horizontal Branch models used in build- 
ing the SSPs (see text) for two different metallicities and several 
masses (0.55 to 1.5 Mq; left to right). Notice the spread in T ef j 
and luminosity, thus the impossibility to treat the HB like a dot. 





3.70 3.57 3.43 3.30 
Log (Teff/K) 



3.70 3.57 3.43 3.30 
Log (Teff/K) 



Figure 9. The evolution along the red giant (dashed line) and 
asymptotic branch (solid line) for a 1 Mq star with solar mctal- 
licity with (right panel) and without mass loss (left panel). In 
the latter case the star evolves to higher luminosities and forms 
a carbon star (diamond) while in the former case its evolution is 
stopped at a lower luminosity. Therefore, a proper modelling of 
mass loss is important, otherwise the synthetic population would 
be too red. 



termined by the star formation process itself. In the next 
two sections we briefly present the main ingredients of the 
synthetic stellar population and of the chemical evolution 
model. The latter is basically the same used to reproduce 
the chemical composition of the Milky Way, since we are 
testing the hypothesis that all LSBs originate in the same 
way as HSBs, apart from being hosted in dark halos with 
larger spin parameter. 



3 SYNTHETIC SPECTRA AND COLOURS 

In order to compute the spectro-photometric evolution of 
LSBs we have used a new set of synthetic stellar population 
models that are an updated version of the previous set of 
Jimenez et al. (1998) models. Our models are based on the 
extensive set of stellar isochrones computed by Jimenez et al. 
(1998) and the set of stellar photospheric models computed 
by Kurucz (1992) and Jimenez et al. (1998). The interior 
models were computed using JMSTAR15 (Jimenez et al. 
1998) which uses the latest OPAL95 radiative opacities for 
temperatures larger than 8000 K, and Alexander's opaci- 
ties (private communication) for those below 8000 K. For 
the stellar photospheres with temperatures below 8000 K 
we have used a set of models computed with an updated 
version of the MARCS code (U. J0rgensen, private commu- 
nication) , where we have included all the relevant molecules 
that contribute to the opacity in the photosphere. For tem- 
peratures larger than 8000 K we have used the set of photo- 
spheric models by Kurucz (1992). Stellar tracks were com- 
puted self-consistently, i.e. the corresponding photospheric 
models were used as boundary conditions for the interior 
models. This procedure has the advantage that the stellar 
spectrum is known at any point on the isochrone and thus 
the interior of the star is computed more accurately than if 
a grey photosphere were used, and therefore we overcome 
the problem of using first a set of interior models computed 
with boundary conditions denned by a grey atmosphere and 
then a separate set of stellar atmospheres, either observed or 
theoretical, that is assigned to the interior isochrone a pos- 
teriori. The problem is most severe if observed spectra are 
used because metallicity, effective temperature and gravity 
are not accurately known and therefore their position in the 
interior isochrone may be completely wrong (in some cases 
the error is larger than 1000 K). A more comprehensive dis- 
cussion and a detailed description of the code can be found 
in Jimenez et al. (1998). 

An important ingredient in our synthetic stellar popu- 
lation models is the novel treatment of all post-main evo- 
lutionary stages that incorporates a realistic distribution of 
mass loss. Thus the horizontal branch is an extended branch 
(see Fig]8|) and not a red clump like in most stellar popula- 
tion models. Also the evolution along the asymptotic giant 
branch is done in a way such that the formation of carbon 
stars is properly predicted (see Fig|j) and therefore the ter- 
mination of the asymptotic giant branch. It is worth noticing 
that in most synthetic stellar population models a constant 
mass loss law is applied resulting in horizontal branches that 
are simply red clumps. This leads to synthetic populations 
that are too red by 0.1 to 0.2 in e.g. B — R (Jimenez et al. 
1998). 

Using the above stellar input we compute synthetic stel- 
lar population models that are consistent with their chemical 
evolution. The first step is to build simple stellar populations 
(SSPs). An SSP is a population of stars formed all at the 
same time and with homogeneous metallicity. The procedure 
to build an SSP is the following: 

(i) A set of stellar tracks of different masses and of the 
metallicity of the SSP is selected from our library. 

(ii) The luminosity, effective temperature and gravity at 
the age of the SSP is extracted from each track. Each set of 
values characterizes a star in the SSP. 
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Figure 10. SSPs for Z = 2Zq. Ages range between 1 X 10 6 and 1.5 X 10 10 yr (from top to bottom). 



(iii) The corre- 
sponding self-consistent photospheric model is assigned to 
each star. 

(iv) The fluxes of all stars are summed up, with weights 
proportional to the stellar initial mass function (IMF). 

SSPs are the building blocks of any arbitrarily compli- 
cated population since the latter can be computed as a sum 
of SSPs, once the star formation rate is provided. In other 
words, the luminosity of a stellar population of age to (since 
the beginning of star formation) can be written as: 



Lx(t ) = 



Lssp,\(Z, t -t)dZdt 



(13) 



where the luminosity of the SSP is: 



L S sp,x(Z,to-t) = / SFR(Z,M,t)l x (Z,M,t -t)dM(U) 

and l\(Z, M, to — i) is the luminosity of a star of mass M, 
metallicity Z and age to — t, Zi and Zf are the initial and final 
metallicities, M m i n and M ma x are the smallest and largest 
stellar mass in the population and SFR(Z, M, i) is the star 
formation rate at the time t when the SSP is formed. 

We have computed a very large number of SSPs, with 
ages between 1 x 10 6 and 1.5 x 10 10 years, and metallicities 



from Z = 0.0002 to Z = 0.1, that is from 0.01 to 5Z e . Fig. |k] 
shows a set of synthetic spectra of SSPs with Z = 2Zq 
and ages between 1 x 10 6 and 1.5 x 10 10 years (from top to 
bottom). 

Our synthetic stellar population models have been ex- 
tensively used in previous works (e.g. Dunlop et al. 1996; 
Spinrad et al. 1997), where a detailed comparison with other 
models in the literature can be found. 



4 CHEMICAL EVOLUTION MODELS 

Most synthetic stellar population codes do not account for 
the chemical time-evolution of the stellar population, and 
use a simple analytical prescription for the star formation 
rate (e.g. the Hubble sequence is often explained like a con- 
volution of SSPs with constant metallicity (solar) and an 
exponentially decaying star formation rate with different e- 
folding times). This approach is not self-consistent, but it 
has been historically used because of the lack of theoretical 
libraries of stellar photospheres of non-solar composition. 
Indeed, the way to proceed is to use chemical evolution mod- 
els where Z(t) and SFR(t) are consistent with each other. 
The model of chemical evolution we use is similar to that 
described in Matteucci & Francois (1989) developed to cal- 
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Figure 11. Time evolution for several colours predicted by our models. Two radii are shown: the nucleus (solid line) and the outermost 
radius of the disk (dotted line). The shaded regions are the observed values for the bluest galaxies in De Blok, Van Der Hulst & Bothun 
(1995) with lcr errors. In all cases the best fit is at late epochs (> 7 Gyr), what shows that LSBs are not young galaxies. 



culate the chemical evolution of the Milky Way. The main 
assumptions of the model are: 

• The disk is represented by several circular concentric 
shells 2 kpc wide and no exchange of gas between them is 
allowed. 

• The disk is formed by infall of primordial material and 
the infall rate is higher in the center than in the outermost 
regions of the disk. In particular, the infall law is expressed 
as: 

5W(r,t) = A(r)X} nf e- t/T(r) (15) 

where r(r) is the timescale for the formation of the disk at 
a radius r. The values of r(r) are chosen in order to fit the 
present time radial distribution of the gas surface density in 
the disk. In analogy with what is required for the disk of 
the Milky Way, we assumed an "inside-out" mechanism of 
formation of such galaxies, implying that r(r) is increasing 
towards larger radii (see Table 2). X\ n f is the abundance of 
element i in the infalling gas and the chemical composition 
is assumed to be primordial. The parameter A(r) is obtained 
by requiring the surface density now to be E(r, t now ), given 
by the disk model described in section 2: 

*M = ril^tlL) (16) 

• The IMF is taken from Scalo (1986) and is assumed to 
be constant in space and time. 

• The star formation rate is assumed to depend on both 
the gas surface density, E gas (r, t) and the total mass surface 
density, E(r, t). In particular, we adopted a law of the type: 

V(r,t) = vZ gaa (r,t) kl Z(r,t) k2 (17) 

with fci = 1.5 and — 0.5. This choice is the best for 



the disk of our galaxy (Tosi 1982; Matteucci 1986). The pa- 
rameter v, which represents the efficiency of star formation 
(the inverse of the time scale for star formation) has been 
assumed to be vi s b = 0.25Gyr~ . 

The effect of stellar winds and outflows is properly ac- 
counted for using the energy output due to supernovae. 
Therefore, at each ring the model accounts for the possi- 
bility of outflow of gas due to the fact that the temperature 
of the interstellar gas is heated enough so the velocity of the 
gas becomes larger than the scape velocity. 

The evolution of several chemical species (H, D, He, C, 
N, O, Ne, Mg, Si, S, Ca, Fe and Zn), as well as the global 
metallicity Z, is followed by taking into account detailed 
nucleosynthesis prescriptions by including the contributions 
to galactic chemical enrichment by stars of all masses. The 
instantaneous recycling approximation was relaxed and the 
stellar lifetimes were taken into account (for the basic equa- 
tions see Matteucci & Francois (1989)). 

In particular, we adopted the Renzini & Voli (1981) 
stellar yields for low and intermediate mass stars (0.8 < 
M/Mq < 8) which produce mostly 4 He, 12 C, 13 C and 14 7V. 
Newer stellar yields for this range of masses are now available 
(Van Den Hoek & Groenewegen 1997) but they do not differ 
significantly from the Renzini & Voli (1981) ones (Matteucci 
et al., in preparation). For the yields of a-elements and Fe 
in massive stars (M > IOMq) we adopted the yields of 
Woosley & Weaver (1995) and for type la supernovae (SNe) 
those of Thielemann, Nomoto & Hashimoto (1996). Type la 
SNe are assumed to be the outcoming of exploding white 
dwarfs in binary systems and their contribution to galac- 
tic chemical evolution is computed in the way described in 
Matteucci & Greggio (1986). They produce mostly Fe and 
some traces of a-elements. For the yields of Zn we assumed 
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Figure 12. Area weighted colours for De Blok, Van Der Hulst & Bothun (1995) sample compared with our model predictions. Also 
plotted the region where colours correspond to ages larger than 7 Gyr. It transpires from the figure that LSBs have ages in excess of 7 
Gyr. Of special interest is the B — I vs B — R plot: it shows a good correlation because it samples the old stellar population and is not 
contaminated by recently formed stars like U. 
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Table 2. Values for the initial surface density (£), t, final metal- 
licity (Z) and final surface density of the gas E ga s for the two 
disc models considered in this work: A = 0.03 (top panel) and 
0.06 (bottom panel). 

the prescriptions given in Matteucci et al. (1993) where this 
element is assumed to originate partly from s-processing in 
massive stars (weak component), partly from r-processing in 
low mass stars (main component) and mostly from explosive 
nucleosynthesis in SNe of type la. 

We computed the evolution of a disk with total final 
mass M d = W 10 M e , Q b = 0.05, VL a = 1 and A = 0.03 and 
0.06 (see Table 2). 



5 THE MODEL GALAXY 

Using the method described in the three previous sections, 
we have computed colours and metallicities for the two 
models described in Table 2. We have assumed a halo for- 
mation redshift z = 0.7, a baryon fraction ft = 0.05, 
a density parameter fio = 1-0 and the Hubble constant 
Ho = 75 kms _1 Mpc _1 . As we said before, we use the sam- 
ple of De Blok, Van Der Hulst & Bothun (1995) to compare 



colors and surface brightness of our models with the obser- 
vations. We also compare our predicted oxygen abundances 
to those measured in McGaugh (1994), and our predicted 
evolution of Zn vs redshift with the data by Pettini et al. 
(1997). 

5.1 Colour Evolution and Ages 

Fig. |ll| shows the time evolution of U-B, B-V, V-I and B-R 
colours for a model with halo mass M = I.elO Mg, halo 
redshift formation z = 0.7, and spin parameter A = 0.06. 
The central radii (solid line) and the outermost radii (dotted 
line) are plotted. 

Also over-plotted as shaded regions are the mean values 
of U-B, B-V, V-I and B-R using F568-1, F571-5, UGC1230, 
F577-V1, UGC128 and UGC628 from De Blok, Van Der 
Hulst & Bothun (1995) for the central region and for the 
outermost radii (the extension of the shaded regions corre- 
sponds to the dispersion among the mean color of the galax- 
ies). We choose this particular set of galaxies because they 
are the bluest in De Blok, Van Der Hulst & Bothun (1995) 
sample and therefore they set the most conservative con- 
straint on the age (lowest age). Two important conclusions 
can be drawn from this plot: the best fit occurs in all cases 
for ages bigger than 7 Gyr, this happens in the nuclear re- 
gion as well as in the outermost radii. Our model provides 
a good fit for both the nucleus and the outer radii simul- 
taneously. It is worth mentioning that the observed colours 
of the nucleus are rather red, unequivocally indicating the 
presence of an old population in the LSBs. This has been 
recently confirmed by Quillen & Pickering (1997) who have 
found a population of old, red stars (B — H > 3.5), with 
colours similar to those of ellipticals, in a sample of LSBs. 

A more interesting constraint on the age can be placed 
by using colour-colour plots. In this case we use the area 
weighted colours for the whole De Blok, Van Der Hulst & 
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Figure 13. Predicted colour profiles for the model considered 
in the text (lines) are compared with the bluest galaxies in De 
Blok, Van Der Hulst & Bothun (1995) sample (dots, see text). 
The agreement is excellent but our profiles arc slightly steeper 
than the observed ones. 

Bothun (1995) sample to compare our predictions. In Fig.|l2| 
we show the evolution of our model (area weighted) for two 
colour plots: U - B vs B - V and B — R vs B — I. We 
have also drawn the lines where the colours correspond to 
ages bigger than 7 Gyr. It transpires from the plots that the 
observations are better fitted if LSBs are older than 7 Gyr. 
The observed points show a better correlation in B — R vs 
B — I than in U — B vs B — V. This is because I and R 
sample much better the old stellar population than U. U is 
affected by recent episodes of star formation, and this is the 
possible reason for the scatter in the U — B vs B — V plot. 
However, there are LSBs with observed U — B > —0.2 and 
B — V > 0.5, which give ages > 7 Gyr, in agreement with 
the B — R vs B — I age determination. 

It is worth noticing that the A = 0.03 model has colours 
that are 0.2 redder than the A = 0.06 model with a similar 
time evolution. 

These more sophisticated models confirm the re- 
sults obtained in our previous paper (Padoan, Jimenez & 
Antonuccio-Delogu 1997), where we were using the most 
simplistic model to show that LSBs have an underlying pop- 
ulation older than 7 Gyr. 

5.2 Colour Profiles, Surface Brightness and 
Spectra 

An interesting test for our model is to see if it can reproduce 
the colour profiles measured by De Blok, Van Der Hulst & 
Bothun (1995). Fig|l| shows our model prediction compared 
with the mean colours for F568-1, F571-5, UGC1230, F577- 
VI, UGC128 and UGC628 from De Blok, Van Der Hulst 
& Bothun (1995) sample with error bars computed as the 
la value of the dispersion among the measured values. The 
agreement is in general rather good. The observed scatter 
in De Blok, Van Der Hulst & Bothun (1995) sample can be 
easily explained by small changes in the spin parameter or 
star formation rate. 

We have also computed the radial surface brightness for 
different times (Fig.fWf). We compare it with the averaged 
values for F571-5, UGC1230, F577-V1 and again the error 



bars correspond to the la value of the dispersion among the 
measured values. Our model is quite successful in fitting the 
observed surface brightness gradients. There is little evolu- 
tion with time in the surface brightness after 7 Gyr, while 
the model is definitely brighter than the observed galaxies, 
for ages smaller than 5 Gyr. 

Fig.[l5] and jlf^ show the predicted spectra for six differ- 
ent ages in the nuclear region (Figjlfj) and in the outermost 
radius (Fig.|l6|). 

5.3 Chemical Abundances 

In the previous sections we have shown that the first stars 
in LSBs formed about 8-9 Gyr ago. One observable which 
can potentially cast doubt on the model is the metal abun- 
dance, for example the Zinc measurements of Pettini et al. 
(1997). In order to ensure that the metal abundance is not 
overpredicted, we present an illustrative computation which 
is designed to maximise the predicted abundance. We as- 
sume a halo collapse at z — 4, corresponding to a look-back 
time of 8 Gyr if Ho = 75 km s _1 Mpc -1 . This extreme 
collapse redshift would correspond to rare high-peak fluc- 
tuations (e.g. Peacock & Heavens 1985), and leads to high 
surface densities, high infall rates and, of course, an early 
build-up of metals. It should therefore be regarded as pro- 
viding a practical upper limit to the elemental abundances 
computed in this section. 

Using our chemical evolution model it is possible to pre- 
dict the evolution with radius and redshift for most chem- 
ical species. In this section we compare the prediction of 
our model with two samples. McGaugh (1994) has observed 
HII regions in a sample of LSBs and derived their oxygen 
abundance. Most of the HII regions are located in the out- 
ermost part of the spiral arms. He finds quite a spread in 
oxygen abundance among LSBs, with a clear peak at about 
[O/H] = -0.9 and dispersion 0.15 (FWHM). His measures 
also show a tentative peak at about [O/H] = —0.4, but as 
is said in McGaugh (1994), it is partly an artifact of the 
method employed to determine the oxygen abundances and 
not a definitive feature of bimodality. In Figjl7| (upper-left 
panel) we show our model predictions for four radii and 
also show the range of oxygen abundances measured by 
McGaugh (1994) as thick horizontal solid lines. It is clear 
that our model is consistent with the McGaugh (1994) mea- 
sures, but does not explain the lowest abundances measured 
by McGaugh (1994). This may be explained by the fact 
that the outermost radius is always close to (or below) the 
Toomre/Kennicutt threshold and it is therefore likely that 
some HII regions will be forming stars in areas of the galaxy 
with metallicity close to primordial values, while the aver- 
age metallicity of the galaxy at that radius (that is what 
we compute) will be higher. In fact, this is well supported 
by the fact that the colour profiles from the sample by De 
Blok, Van Der Hulst & Bothun (1995) always show redder 
colours in the outermost radius than the ones corresponding 
to a galaxy with [O/H] as low as —1.2. For the outermost 
radius our model predicts [O/H] = —0.8 in excellent agree- 
ment with the peak found by McGaugh (1994). It should 
be also noted that none of the innermost radii in our model 
reproduces McGaugh (1994) measures, in good agreement 
with the well known fact that HII regions are usually in the 
outermost parts of the hosting galaxy. 
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Figure 14. The surface brightness profiles for our model at different ages are compared to the corresponding averaged LSB from De 
Blok, Van Der Hulst & Bothun (1995) sample (see text). It is clear that our model gives a good fit both for the profile and the absolute 
value of the surface brightness in both bands. It also transpires from the figure that ages as low as 4 Gyr do not provide a good fit, 
reinforcing the old age argument presented in this work. 



Another important prediction of our model is the evo- 
lution of Zn. The right-upper panel of Fig.^8| shows the pre- 
dicted evolution of Zn for four different radii in our model. 
We also show as diamonds the data for Damped Lyman-a 
(DLA) systems observed by Pettini et al. (1997), including 
some upper limits. While it is difficult to establish a corre- 
lation of Zn vs redshift, it transpires from the figure that 
the predicted Zn abundance for LSBs at different radii is 
consistent with the Pettini et al. (1997) data. In fact the 
spread found at every redshift is consistent with the spread 
in Zn abundance predicted for LSBs. It should be noted 
though, that this does not mean that DLAs are LSBs, this 
only means that due to their old age and their presence at 
high-z, some DLAs could be caused by LSBs in the line of 



sight. Indeed, at least one DLA has been identified as a LSB 
galaxy (Steidel et al. 1994). 

Notice that the results derived in this section are robust 
for z — regardless of the formation redshift adopted for 
LSBs since the evolution of the chemical elements with time 
is very weak after the first 4 Gyr. 



6 DISCUSSION AND CONCLUSIONS 
6.1 The Nature of LSBs 

In this work we have shown that the photometric proper- 
ties of the bluest galaxies in the sample by De Blok, Van 
Der Hulst & Bothun (1995) can be reproduced assuming 
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Figure 15. Predicted spectra for the nucleus of a typical LSB. 



that they form and evolve as normal disk galaxies, with rel- 
atively high spin parameter. That the spin parameter alone 
could explain the low surface brightness was already shown 
by Dalcanton, Spergel & Summers (1997), but they assumed 
a given M/L ratio. Here we have strengthened the point 
by also explaining the colours, the colour gradients, and the 
chemical abundances. We have also provided synthetic spec- 
tra that could be compared with the observations, in the 
effort to prove that LSB disk galaxies are indeed normal 
galaxies with high spin parameter. The stellar populations 
of LSBs are rather evolved, especially in the central part of 
their disks. Even the bluest galaxies in the LSB sample by De 
Blok, Van Der Hulst & Bothun (1995) are older than 7 Gyr, 
and typically star formation in these galaxies starts about 
9 Gyr ago. This confirms the results of Padoan, Jimenez & 
Antonuccio-Delogu (1997), where the galaxies in the same 
sample were found to be at least 7 Gyr old. 

LSB and HSB disks are therefore hosted in similar dark 
matter halos that differ only for their spin parameter, and 
they evolve in a similar way, apart from the fact that LSB 
disks are less concentrated than HSB disks. This means that 



the present day (or low redshift) abundance of LSBs relative 
to HSBs should be about the same at any redshift. More- 
over, the conclusions about the epoch of star formation and 
halo formation (see the next section), and about the colour 
evolution must be valid for disk galaxies in general. In par- 
ticular, Fig. jll] shows that the colour evolution of LSBs is 
very strong. The U-B, B-V and V-I colours of the central 
part of the disk at high redshift are about 1 mag bluer than 
they are at low redshift, and the B-R colour even 1.5 mag 
bluer. This must be a property of disks in general, and not 
only of LSB disks. 



6.2 Star Formation and Galaxy Formation 

We have shown in section 2.1 that the size of LSB disk galax- 
ies is such that their halos are formed in the redshift range 
0.1 < z < 1.2. We have also been able to estimate the for- 
mation redshift of the halo of each galaxy in the sample, 
and we have found a mean value z = 0.7. In an Einstein- 
de-Sitter universe with Ho = 75 kms _1 Mpc _1 , a redshift 
z = 0.7 corresponds to a look-back time of 6.8 Gyr. On 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy Formation and Evolution: Low Surface Brightness Galaxies 15 




Figure 16. Predicted spectra for the outermost radius of a typical LSB. 



the other hand, the photometric properties of LSB galax- 
ies, in particular their colours, are such that star formation 
must have started approximately 9 Gyr. Therefore, the star 
formation starts about 2 Gyr before the galactic halos are 
formed. These numbers are likely to be similar for HSB disk 
galaxies too, since they have colours and sizes very similar 
to LSBs, also suggesting a rather low formation redshift for 
their halos. 

The discrepancy between the epoch when star forma- 
tion starts in a galaxy and the epoch when the halo forms 
means that the process of star formation starts before the 
galactic halo has assembled, as predicted for example by 
bottom-up galaxy formation scenarios, characterized by 
power spectra with more power on sub-galactic scales than 
on galactic ones. 

In Fig. ^ we show the evolution with redshift of the 
SFR. The upper panel shows the fraction of gas with time 
that is converted into stars at every redshift for the nucleus 
and the disk. The continuous thick line shows the average 
of the nucleus and the disk as a representative of the whole 
LSB. It is worth noticing that most of the gas in the nu- 



cleus is converted into stars between 2 < z < 4, while most 
of the star formation in the disk takes place at z < 2. The 
lower panel shows the global SFR in the Universe for LSBs 
assuming that they have the same comoving number density 
as HSBs and that the typical mass for a L* LSB is a fac- 
tor 5 smaller than the typical mass of a normal L, galaxy. 
The data points are taken from Madau (1997) and trace the 
global SFR in the Universe for HSBs. LSBs lie below the 
observed points because they have typical masses smaller 
than HSBs. The important point to notice here is that the 
global SFR for LSBs has a flatter evolution than the ob- 
served points. 

6.3 Discussion 

In order to explain the rotational properties of LSBs, we 
require halos with velocity dispersions w 100 km s" 1 , or 
masses w 10 11 M©. This is also consistent with the required 
gas masses and a baryon fraction of a few percent. To obtain 
the right sizes of LSBs, we additionally require a relatively 
late halo collapse, z ~ 0.7. Clearly this may present a prob- 
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Figure 17. The redshift evolution for some elements is plotted. The top-left panel shows the evolution of [O/H] and is compared with 
McGaugh (1994) measures. Our model reproduces perfectly the peak in [O/H] among LSBs found by McGaugh (1994) at [O/H] = —0.9. 
The top-right panel shows the analogous to the previous one but for [Zn/H] compared with Pcttini ct al. (1997) measures. It transpires 
from the figure that the predicted [Zn/H] for LSBs is consistent with the spread found for DLAs. 



lem for galaxy formation models which have a lot of power 
on the appropriate scale. An obvious example is COBE- 
normalised CDM, with shape parameter F = 0.5 (Frenk 
& White 1991), for which the characteristic collapse red- 
shift for 100 km s _1 halos is around z = 6. Reducing the 
amplitude to that of standard CDM (bias parameter 2.5) 
alleviates the problem to a certain extent, but the collapse 
redshift is still z = 3. It is more promising to consider mixed 
dark matter models, for which the power spectrum is in any 
case a better match to the observed galaxy spectrum. For 
a neutrino fraction of 30%, the comoving number density 
of halos with mass exceeding 1.5 x 10 11 M reaches 0.01 
h 3 Mpc~ 3 , comparable to the density of large galaxies, at 
a redshift of about 1.7 Klypin et al. (1995). This is still a 
little high, but it should be remembered that the collapse 
of halos of given mass does take place over a reasonably 
wide range of redshifts, and some fraction of halos would 
form after z = 0.7. In the model considered by Klypin et al. 
(1995), about a quarter of halos exceeding this mass form 
after z = 0.7. Analytically, there is expected to be a weak 
correlation between high spin and late collapse Heavens & 



Peacock (1988), and numerical simulations by Ueda et al. 
(1994) show an anticorrelation between spin and density. 

We postulate therefore that the LSBs form in halos with 
high spin which form relatively late, in the context of a hi- 
erarchical model. The indications are that a model with rel- 
atively little small-scale power, such as mixed dark matter, 
would do best here, but there are enough uncertainties in 
the collapse redshift, both in the simplifications of the halo 
model, and in physical processes which might delay collapse 
(e.g. Babul & Rees 1992) that this conclusion is not strong. 

In summary, the most important results of this work 

are: 

(i) Observed colour profiles, chemical abundances and 
surface brightness profiles for LSBs are well fitted if they 
are assumed to have a spin parameter for its halo higher 
than HSBs. 

(ii) LSBs are not young objects, as often claimed in the 
literature, since their colours are well fitted by old (> 7 Gyr) 
stellar populations. 

(iii) There is discrepancy between the photometric age of 
the galaxies and the age of formation of their halo, which in- 
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Figure 18. Top panel shows the gas fraction with time in our 
LSB model that is converted into stars vs redshift for the disk 
and the nucleus. The thick line corresponds to the average of the 
nucleus and the disk and should be a fair representative of the 
whole LSB. The bottom panel shows the predicted global SFR in 
the Universe due to LSBs assuming that the mass for a L* LSB 
is 5 smaller than the mass for a L* HSB and that the comoving 
number density for LSBs is the same as for HSBs (see text). The 
observed points are taken from Madau (1997). Notice that the 
evolution with redshift is flatter than the one measured by Madau 
(1997). 



dicates that star formation can start about 2 Gyr before the 
halo is formed. This is perfectly acceptable in hierarchical 
models of galaxy formation, and, as discussed in the text, 
the discrepancy can be reduced or removed if the cosmo- 
logical model we have assumed is incorrect. If an open (A) 
cosmology had been adopted, the above discrepancy would 
have been completely removed. 
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